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Mixed  Basis  Band  Structure  Interpolation  Scheme  Applied  to 
the  Fluorite  Structure  Compound*  NiSi2>  AuA12,  AuG*2  .  and  Auln2 

Sehun  kirn,  Jeffrey  G.  Nelson,  and  R.  Stanley  Williams 

I 

Department  of  Chaistry  and  Biochemistry 
University  of  California  Los  Angeles 
Los  Angeles,  California  90024 


A  mixed  basis  band  structure  interpolation  scheme  for  fee  d-band 

t  ■ 

petals  has  been  extended  to  include  fluorite  structure  (CeF^j  compounds  by 
incorporating  more  plane  waves  in  the  basis  set.  Since  the  fluorite  and 
icc  structures  belong  to  the  same  space  group,  the  interpol etion  scheme 
originally  developed  for  fee  d-band  metals  is  also  capable  of  pnerating 
iiucrite  bend  structures.  Ihe  interpolation  parameters  for  NiSi^  AuA12. 

AuGa^  an(j  Auln2  have  been  determined  by  fitting  nonrela  tiv  istic  first 
principles  calculations  using  a  nonlinear  least  squares  procedure.  Good 
agreement  with  the  first  principles  results  is  obtained  up  to  about  5  eV 
above  the  Fermi  level  for  a  basis  set  containing  39  plane  waves  and  5 
d-functions.  The  parameters  for  the  intermetal  1  ic  compounds  containing  Au 
were  then  adjusted  to  include  the  effects  of  spin-orbit  splitting  in  the 
d-bands  and  to  improve  the  agreement  of  the  calculated  de nsi ty-of-states 
with  the  results  of  photoelectron  spectra.  The  adjusted  d-bands  of  AuAl2> 

Aj'ia2  and  Au  In2  differ  considerably  frem  those  calculated  by  first 

principles.  \  I  \  v  ;  •  v  A-  r  3L  ,  _  ,  v 

,  ,  ,  '  v"'!  '  *  K. V 


I.  Introduction 


Mired  basis  band  structure  interpolation  scbeacs  have  proven  to  be 
extremely  valuable  computational  tools  for  the  study  of  the  electronic 
structure  of  d-band  metals.  The  first  routines  were  developed  by  Hodges, 
hhrenreich,  and  Lang  and  Mueller  for  simple  fee  d-band  metals.  Smith  and 
coworkers'*  ''  have  continued  the  develoiment  of  a  fee  d-band  interpolation 
scheme  over  the  past  decade  and  have  applied  it  as  an  aid  in  interpreting 
nkotoem ission  and  UV  reflectance  data.  The  advantages  in  using  such  a 
second  principles  technique  are  1.)  it  is  simple  and  ineapensive  to 
generate  E  vs.  k  and  densities-of-atates  plots,  2.)  the  parameters  in  an 
interpolation  scheme  may  be  easily  adjusted  to  improve  the  agreement 
between  the  calculated  energy  bands  and  experimental  observations,  and  3.) 
the  parameters  may  be  reported  in  the  literature  and  used  by  other 
investigators  to  generate  the  identical  band  structure,  which  makes  the 
results  of  a  band  structure  de  terming  tion  very  portable.  Furthermore, 
Smith  and  coworkers  have  shown  that  transi tion  ma tr ix  elements  calculated 
using  the  parameterized  Hamiltonian  agree  reasonably  well  with  observed 
spectral  intensities  in  reflectance  and  Angle  Resolved  Photoemission 
Spectroscopy  (ARPES).^'®  By  fitting  a  band  structure  interpolation  scheme 
to  experimental  data,  it  should  be  possible  to  experimentally  determine  the 
electronic  structure  of  solids  to  a  high  degree  of  accuracy. 

In  this  paper,  the  interpolation  scheme  of  Smith,  et  al .  **  is 

extended  further  to  more  complex  systems,  i.e.,  the  fluorite  structure 
compounds  Afi^ ,  which  have  larger  and  more  complex  unit  cells  than  the  fee 
d-band  metals.  The  fluorite  structure  is  composed  of  one  A  and  two  B  fee 
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2 


sub  lattices,  with  one  B  sublattice  translated  by  one  fourth  the  body 
diagonal  of  the  cubic  unit  cell  along  both  the  [111]  and  [Ill]  directions 


%  with  respect  to  the  A  sublattice.  Several  intermetal  1  ic  compounds 

involving  Au  and  Pt  (the  A  elements  of  AB^ )  and  the  technological  iy 
important  silicides  NiSi^  and  CoS»2  have  the  fluorite  structure.  In  these 
cases,  the  A  atoms  have  valence  d-orbitals  that  are  very  important  in 
determining  the  electronic  structure  of  the  compounds,  but  the  d 
contribution  of  the  B  atoms  is  probably  extremely  small.  Since  the  A  atoms 
reside  on  a  fee  sublattice  with  larger  dimensions  than  in  the  corresponding 
elements,  the  d-d  interactions  in  the  compounds  should  be  much  smaller  than 
ir.  the  elements.  Thus,  it  should  be  quite  interesting  to  compare  the 
d-band  regions  of  the  various  intermetal  1 ic  compounds  with  each  other  and 
wiiii  the  corresponding  fee  d-band  metals.  However,  to  date  very  few  band 
structure  calculations  have  been  performed  on  the  d-band  fluorite 
compounds . 


At  first,  the  application  of  an  interpolation  scheme  to  the  fluorite 
|  compounds  may  seem  complicated,  but  since  the  fluorite  structure  has  the 

same  space  group  symmetry  as  fee  crystals,  the  mixed  interpolation  scheme 

developed  for  fee  d-band  metals  is  also  quite  satisfactory  for  the  fluorite 

L 

structure.  An  interpolation  scheme  has  been  developed  similar  to  that  of 

:  t  h ,  et  al.,3  ®  except  that  more  plane  waves  were  included  in  the  basis  _ _ 

ct  to  enhance  the  convergence  in  the  fitting  procedure.  Single-zeta 
-ter  orbitals  were  used  to  approximate  the  radial  part  of  the  d-wave 
motions  in  calculating  d-wave  overlap  integrals,  but  this  approximation 
?>  no  better  than  Smith's  use  of  spherical  Bessel  functions  for  the 

Using  a  nonlinear  least  squares  procedure, ^ 
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interpolated  Bands  may  be  fit  to  the  energy  eigenvalues  at  high  symmetry 
points  in  the  Brillouin  Zone  ( BZ)  of  first  principles  calculations  to 
determine  the  nineteen  parameters  used  in  the  interpolation  scheme.  The 
resulting  band  structures  may  be  optimized  by  adjusting  some  of  the  fitting 
parameters  to  agree  with  experimental  observations,  and  then  compared  with 
one  another  to  analyze  the  effects  of  atcmic  structure  and  composition  on 
electronic  structure. 


The  main  purpose  of  the  research  reported  here  was  to  determine  how 

well  a  simple  interpolation  scheme  could  compute  the  band  structure  of  the 

d-band  lluorite  compounds.  In  Sec.  II,  the  details  of  the  present 

interpolation  scheme  are  summarized.  These  include  a  description  of  the 

3-b 

mixed  interpolation  scheme  of  Smith,  et  al.,  the  modification  of  the 

basis  set,  the  single-zeta  type  overlap  integrals,  and  a  snamary  of  the 

7 

nonlinear  least  squares  procedure.  In  Sec.  Ill,  the  results  of  the 

application  of  the  present  scheme  to  fit  nonrel  a  tiv  ist  i  c  first  principles 

band  structure  calculations  for  NiSij,  AuAl^,  AuCa2>  and  Auln^  are 
presented.  Finally,  Sec.  IV  contains  a  discussion  of  these  results,  an 
extension  of  the  calculations  to  include  spin-orbit  effects  and  improve 
agreement  with  experimental  results  in  the  Au  Sd  bands  of  AuAl^,  AuGa^  and 
Auln^,  and  a  further  comparison  of  the  results  to  experimental  data  and 
band  structures  of  the  fee  metals,  Ni  and  Au. 


r  v  •'$* 


interpolation  Scheme  and  Fitting  Procedure 


This  work  basically  adopts  the  aixed  interpolation  scheme  developed  by 
Smith,  et  al.3-6  According  to  this  scheme,  the  nonrel ativ istic  Hamiltonian 
ior  tcc  d-band  metals  is 


H  H  A 

cc  cd 

*aa 

dc  dd 


(1) 


*i'sre  ^cc>  H  ^  e°d  represent  the  plane  wave,  the  hybridisation  and  the 

d-orbi tal  tight  binding  blocks,  respectively. 

All  of  these  blocks  retain  the  saae  parameterized  forms  as  in  the 
i cterpol  a ti  on  scheme  of  Saith,  et  al.3  ®  However,  the  plane  wave  block  of 

has  been  extended  to  a  39  x  39  rather  than  16  x  16  sub-aatrix  in  order 

c  c 

to  reproduce  the  valence  band  structure  of  the  fluorite  coapounds.  In 
terms  of  the  reduced  wave  vector  k,  the  39  plane  waves  included  in  the 
basis  set  are  those  with  wave  vectors  k.  *  k  ♦  G.,  where  k  is  constrained 
tc  be  within  the  first  BZ  and  includes  the  reciprocal  lattice  vector 
'0,0,0),  all  eight  vectors  of  the  type  ( 2n / a ) (+1 ,+l ,+l ) ,  all  six  of  the 
type  (2n/a) (*2,0,0) ,  all  twelve  of  the  type  (2n/a ) (+2 ,+2 ,0)  and  all  twelve 
.:(  the  type  (2n/a)(-3»*l,+l). 

The  Hamiltonian  elements  of  the  plane  wave  block,  H  are  expressed 

-  r,  the  same  form  as  in  Smith's  scheme. 
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The  first  tern  is  the  free  electron  energy,  in  which  a  =  (-h^/2m)  (n/4a)^  and 


a  is  the  lattice  constant.  The  V  are  local  pseudopotential  tens,  of 
which  the  lowest  six,  V000,  Vlll,  V200,  V220,  V311  and  V222  are  assumed  to 
be  nonzero.  The  paraaeter  S  is  a  constant  in  the  or thogonal  ization  tern, 
j2  is  the  usual  spherical  Bessel  function  of  order  2,  and  is  the 

Legendre  polynoaial  of  order  of  2.  The  S  ten  aay  be  regarded  as  a 
nonlocal  pseudopotenti al  and  the  paraaeter  R  ps  a  Buff  in- tin  radius;  both 
are  treated  as  disposable  paraaeters  for  this  scheae. 


The  fora  of  the  local  pseudopotential  in  the  plane  wave  portion  of  the 
parameterized  Hamiltonian  for  the  fluorite  structure  aay  be  expected  to 
differ  from  that  for  fee  aetals.  The  total  local  pseudopotential  of  AB2 
can  be  written  as  follows: 


Vg(AB2)  =  Vg<A)  ♦  Vg(B)cos(G*x).  (3) 

where  t= ( a / 4 ) ( 1 , 1 ,1 )  and  x  is  a  vector  connecting  an  A  ataa  to  one  of  the  B 
atoms  in  the  primitive  unit  cell  (the  other  B  atom  is  at  -x)  of  the 
fluorite  lattice  structure.  The  VQ(A)  and  VC(B)  pseudopotential  tens  are 
in  principle  dependent  only  on  the  identity  of  the  A  and  B  atoms, 
respectively.  For  the  G  vectors  of  the  reciprocal  lattice  of  the  fluorite 
structure,  the  values  of  cos(G.*  x)  are  1  iaited  to-1,  0,  or  1  and  are 
determined  only  by  the  Bagnitude  of  G.,  not  its  direction.  Although 
equation  (3)  is  expressed  as  a  sub  of  the  pseudopotential  coefficients  for 
two  different  atoms,  the  total  local  pseudo po tenti al  has  exactly  the  same 
functional  lorm  as  the  pseudopotential  parameter' s  in  Saith*  s  fee  d-band 
sth  eme. 
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Thus,  the  values  of  >G  for  ABj  *re  ell  constants,  with 


VG  (A)  "  VG  (B)  for  |G.|  -  2,  (71 

~  i  ~i 


VG  (AB2>  =  VG  (A) 
~  i  ~  i 


for  |Ga  (4) 


VG  (A)  ♦  VG  (B)  for  |G i |  =  8,  4, 

~  i  ~i 


where  I G  ^ |  is  in  units  of  2n/a.  As  long  as  the  d-orbital  contribution  of 
the  B  atoms  to  the  band  structure  is  negligible,  the  interpolation 
Hamiltonian  of  the  fluorite  structure  is  identical  to  that  of  the  fee 
structure. 


The  Hamiltonian  elements  of  the  hybridization  block,  H  ,,  are 

cd 

expressed  in  the  following  form: 


“u  ■  Bt.«viR>yy- 


where  B ?  and  are  constants  in  the  hybridisation  tanas  that  are  treated 

as  disposable  parameters  for  orbitals  with  and  e^  symmetry, 

respectively  The  functions  Y^(k^)  represent  the  set  of  real  spherical 

harmonics.  The  spherical  Bessel  functions  in  Eqs.  (2)  and  (5)  are  overlap 

integrals  involving  the  radial  parts  of  the  d-orbital s  and  plane  waves  in 

2 

which  the  function  r  i*  approximated  by  a  delta  function.  A 

possible  improvement  to  the  overlap  integral  might  be  to  approximate  K  ,(r) 

DO 

with  Slater  wave  functions.  The  resulting  overlap  integrals  nsed  in  this 

work  are  listed  in  Appendix  A.  In  the  fitting  procedure,  the  initial  value 

t i  the  orbital  exponent,  seta,  was  taken  from  the  tables  of  Clementi  and 

8 

Koetti  for  the  appropriate  atomic  d-wave  function,  after  which  seta  was 


■i 


%  '  ?  j^J^;  ^r  t.  ^ 


treated  as  a  disposable  parameter  in  the  fitting  procedure. 


The  matrix  elements  of  the  block  may  be  expreased  in  the 

a 

tight-binding  forms  given  by  Slater  and  Koster  or  Fletcher  and 
Wohl far th. 1 0  The  nearest  neighbor  three-center  form  of  Fletcher11  was  used 
in  the  present  scheme.  The  Fletcher  notation11  involves  eight  parameters  : 

EO ,  »  Aj ,  A^ .  Aj ,  A^ ,  Aj  and  A^ .  The  A  parameters  specify  the  dispersion 

of  the  bands,  while  ED  and  EO  +  A  are  the  mean  energies  of  the  t^  and  e 
subbands,  respectively. 


For  relativistic  band  structure  calculations,  the  effect  of  spin  orbit 

12 

coupling  on  the  d-bands  may  be  included  in  the  manner  of  Friedel,  et  al . 
The  total  Hamiltonian  has  the  following  form: 


"m-w 


where  M  and  N  are  given  by  Friedel,  et  al •  The  spin-orbit  splitting 
parameter  {  may  be  determined  from  ARPES  measurements11  or  be  taken  from 
atomic  values. 


The  parameters  of  the  mixed  interpolation  scheme  for  a  particular 

7 

system  were  determined  by  applying  a  nonlinear  least  squares  procedure  to 

tit  the  energy  eigenvalues  of  first  principles  calculations.1*-1**  In  this 

fitting  procedure,  each  energy  eigenvalue  was  given  equal  weight  regardless 

3  17 

of  its  degeneracy.  The  parameters  of  Ni  and  Au  given  by  Smith,  et  al . ,  ' 


mere  used  as  an  initial  guess  of  the  parameters  for  the  fluorite  structure 


compounds,  except  that  the  a  and  R  parameters  were  changed  to  reflect  the 

different  lattice  constants  of  the  intermetal 1  ic  compounds,  where  a  is 
2 

proportional  to  1/a  and  R  is  expressed  as  a  fraction  of  a.  Using  this 
procedure,  the  parameters  usually  converged  to  values  which  yielded  a  root 
ir.ean  squre  (rms)  error  in  the  energy  bands  of  the  fluorite  structure 
compounds  of  about  0.15  eV  after  3-4  iterations. 

After  the  fitting  parameters  for  a  band  structure  were  determined,  the 

total  densi  ty-of-states  (TDOS)  of  the  system  was  obtained  by  summing  a  set 

of  Gaussian  peaks  with  a  0.3  eV  FW HR  that  were  centered  at  the  energy 

18 

eigenvalues  calculated  at  60  special  points  in  the  BZ.  These  special 

points  were  chosen  to  efficiently  calculate  averages  over  the  BZ  of 

18 

periodic  functions  of  the  wave  vector.  The  TDOS  plots  may  then  be 
ccrpared  to  x~ray  photoemission  (IPS)  data  or  other  probes  which  sample  an 
average  over  the  valence  or  conduction  bands. 

All  the  calculations  were  performed  using  a  VAX  11/780  computer.  The 

CPU  time  of  the  VAX  required  to  produce  the  eigenvalues  at  a  single  k-point 

for  nonxelativistic  and  relativistic  calculations  was  3.7  and  33  sec., 

respectively,  for  the  cases  in  which  the  B  blocks  contained  39  plane 

c  c 


waves. 
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The  number  of  first  principles  band  structures  for  d~band  fluorite 
compounds  is  limited  to  nonr  el  a  ti  v  ist  i  c  Augmented  Plane  Wave  (APW) 
calculations  for  Au  interne  tal  1  ic  compounds  with  AJ ,  Ga  and  In1*  and 
N  i  S  . 1  ^  ‘ 1  **  The  energy  bands  generated  by  the  interpolation  scheme  and  the 
TDOS  for  NiSi^,  AuAl^ ,  AuGa^  and  Auln^  are  illustrated  in  Figs. 1-4, 
respectively.  The  fitting  procedure  utilized  33-39  energy  eigenvalues 
taken  from  the  AiV  results1*  1  at  4-5  symmetry  points. 

The  main  difference  between  NiSi^  and  the  Au  alloys  is  that  in  NiSi^, 
the  d-bands  mix  strongly  with  the  plane  wave  states,  whereas  the  d-bands  in 
the  Au  alloys  essentially  reside  in  a  band  gap  of  the  plane  wave  states. 
One  especially  noteworthy  feature  associated  with  AuGa^  is  the  flat  band 
with  A2  symmetry  calculated  by  Switendick  and  Narath*  s  AW  procedure1*  and 
reproduced  by  the  interpolation  scheme  in  Fig.  3.  In  Figs.  5  and  6,  the 
effects  of  including  spin-orbit  coupling  in  the  d-bands  of  AuGa^  and  Aul^ 
as  well  as  adjusting  the  d-bands  to  agree  with  ARPES  data1^  by  varying  the 
BO  and  A  parameters  are  shown. 

In  Fig.  7,  the  TDOS  for  AuAl^,  AuGa^  and  Auln^  are  compared  with  the 

1 9 

corresponding  x-ray  photoelectron  spectroscopy  (XPS)  data.  In  this 
comparison,  the  raw  experimental  data  with  no  background  subtraction  are 
displayed.  For  each  compound,  two  TDOS  curves  calculated  using  the 
interpolation  scheme  are  shown:  one  for  the  fit  to  first  principle’s 
results  and  the  other  for  the  empirically  adjusted  d-bands  (for  AuAl^ ,  the 
AuOa^  spin-orbit  coupling  parameter  was  used  and  the  BO  and  A  parametrs 
were  adjusted  to  enhance  the  agreement  with  the  XPS  data  of  Ref.  (19] ).  A 

10 


us  s',  an  broadening  u  1  1.0  eV  was  applied  to  all  these  TDOS  curves  to 


>  g  *  1  y  simulate  the  total  experimental  resolution  of  the  spectra. 

The  parameters  determined  from  the  fitting  procedure  for  the  fluorite 
-.pounds  are  listed  in  Table  I,  along  with  the  standard  deviation  for  each 
rsr-.eter  determined  by  the  fitting  procedure.  For  the  purpose  of 

ririson,  this  table  also  lists  the  parameters  of  Ni  and  Au  calculated  by 

scheme  for  39  plane  waves  in  the  U  block.  Some  of  the  A  parameters 

cc 

the  fluorite  compounds  have  negative  values,  unlike  those  of  the  fee 
.i.-Lsis,  but  in  most  of  these  cases  the  uncertainty  in  determining  the 
stive  values  is  as  large  as  those  A  parameters  (i.e.,  they  are 
r.  tively  zero).  Both  the  B  hybridization  and  S  or thogonal  ization 
u iters  of  the  intermetal  1 ic  compounds  are  smaller  than  those  of  the  fee 
c .  s.  The  pseudopotential  coefficients  V200,  V220  and  V311  of  the 

rounds  are  larger  and  V222  is  smaller  than  those  of  the  corresponding 
.metal  s,  which  is  the  result  of  the  contribution  to  the  pseudopotenti  al 
the  Si  or  Group  III  atoms  in  the  compounds.  As  expected  from  Sec.  II, 
values  of  Vlll  for  NiSi^  and  Ni  are  nearly  identical,  since  the  form 
ter  of  the  pseudopotential  of  the  B  atoms  in  AB^  is  zero.  However,  for 
ether  cases  involving  NiSi^  and  the  Au  intense tal  1  ics,  the  values  of 
1  and  V311  differ  substantially. 

The  entries  at  the  bottom  of  Table  I  indicate  the  rms  errors  and  the 
■rrge  deviations  for  each  of  the  fits.  The  rms  errors  in  the  fits  to  the 
and  Au  energy  bands  are  larger  than  those  for  NiSi^  and  Au  alloys 
ause  more  AF*  energy  eigenvalues2®’21  (51  for  Ni  and  90  for  Au)  were 
ed  in  the  fitting  procedure.  The  rms  errors  in  the  fits  to  the  energy 


K 
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bands  ol  the  fee  metals  using  the  present  39  OW  scheme  were  less  than 
those  using  Smith’ s  original  16  OF*  scheme.  Thus,  increasing  the  ntmber  of 
ilsne  waves  in  the  basis  set  enables  a  substantial  improvement  in  the  fit 
to  a  band  structure  without  increasing  the  number  of  fitting  parameters. 
Finally.  the  d-orbital  contributions  to  the  wave  function  of  several  bands 
at  the  [  and  X  points  of  the  BZ  obtained  from  the  eigenvectors  of  the 
parametrized  Hamiltonian  are  listed  in  Table  II.  along  with  those  of 
Bylander,  et  al.*^,  for  NiSi,. 


'M 
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Y .  Discussion 

As  mentioned  in  Sec.  II,  the  purpose  of  increasing  the  mnber  of  plane 

aves  in  the  basis  set  of  the  interpolation  scheme  from  the  16  used  by 
2 

r.ith  to  39  is  to  recover  the  flat  band  in  the  AuGa^  valence  band 

r.d  improve  the  fit  to  the  higher  energy  conduction  bands.  The  rms  error 

2 

cteined  after  convergence  of  the  original  16  OPV  scheme  of  Smith  to  36 
-crgy  levels  of  AuGa^**  at  5  symmetry  points  was  0.54  eV,  and  the  flat 
i  nc  ( f  j  » —  X  3  )  was  not  well  reproduced.  When  the  39  plane  wave  scheme  was 
;cc  to  fit  the  same  set  of  energy  eigenvalues,  the  rms  error  quickly 
merged  to  0.115  eV  within  4  iterations.  However,  even  with  39  plane 
ves,  the  second  band  with  and  the  first  with  symmetry  above  the 
-mi  level  Ep  deviate  1-2  eV  frcm  the  results  of  the  first  principles 
. .calations  for  all  four  compounds.  In  order  to  improve  the  fit  to  these 
i  other  conduction  bands,  V400  and  V420  pseudopotential  terms  were 
luded  in  the  Hamiltonian,  but  they  had  a  negligible  effect  on  the 
iity  of  the  fits.  Improving  the  fit  to  these  conduction  bands  may 
ui re  more  plane  waves  in  the  basis,  more  terms  in  the  Hamiltonian  (for 
.tance  to  explicitly  account  for  the  presence  of  the  d-orbitals  of  Ga  or 
in  AuGa^  and  Auln2).  or  both  extentions.  The  above  mentioned  an^  fj 
Juction  band  eigenvalues  were  omitted  in  the  final  fitting  procedure  for 
lour  intermetal  1 ics,  since  including  them  in  the  fit  increased  the  rms 
c-  greatly,  resulted  in  unreasonable  large  values  of  the  standard 


'  ’alien  of  several  of  the  fitting  parameters,  and  distorted  the  valence 


Another  possible  improvement  would  be  in  the  plane  wave  -  Sd  overlap 
integrals,  which  one  would  not  expect  to  be  well  approximated  by  simple 
spherical  Bessel  functions.  An  attempt  was  made  to  use  single-zeta  Slater 
functions  to  approximate  the  radial  portion  of  the  Sd  orbitals,  as 
described  in  section  II.  For  NiSi^,  the  rms  error  in  the  fit  to  the  AW 
results  using  the  39  OW  scheme  after  3  iterations  improved  slightly  when 
compared  to  the  spherical  Bessel  function  approximation  to  0.133  from  0.146 
eV  without  including  the  and  f  j  conduction  band  points,  but  when  these 
points  were  included,  the  same  difficulty  of  convergence  was  encountered, 
it  was  found  that  there  were  no  negative  A  values  using  this  single  zeta 
approximation. 

In  the  case  of  AuGa^,  the  single  seta  wave  functions  were  not 
effective  in  achieving  a  fit,  even  if  the  J and  conduction  band  points 

*■  -v.  “ 

were  not  included.  The  rms  error  for  this  trial  was  0.335  eV  after  the 
third  iteration,  which  is  about  three  times  larger  than  that  of  the  Bessel 
function  approximation.  This  is  perhaps  contrary  to  initial  expectations 
in  view  of  the  severe  approximation  of  using  the  spherical  Bessel  functions 
for  overlap  integrals.  In  fact,  no  Improvement  in  the  accuracy  of  the  fit 
occuxed  when  the  Slater  wave  function  approximations  was  used.  This  may  be 
because  both  overlap  approximations  actually  have  very  similar  k  dependence 
for  small  k,  but  both  may  still  be  poor  approximations  to  describe  the 
exact  overlap  integrals. 

It  is  useful  to  compare  the  energy  band  results  for  the  fluorite 
compounds  among  each  other  and  with  those  for  the  fee  d-metals.  A  common 
measure  of  the  d-band  width,  EtXjJ-E^),  is  4.8  eV  for  NiSi^,  which  is 


14 


larger  than  the  4.4  eV  bandwidth  of  paramagnetic  fee  Ni  even  though  the 
larger  spacings  between  Ni  atns  in  the  intermetal  1  ic  compound  might  be 
expected  to  produce  a  d-band  narrowing.  In  fact,  examination  of  Table  II 
shows  that  a  significant  amount  of  d-orbital  character  it  present  in  bands 
separated  by  almost  10  eV.  Table  II  also  shows  that  the  d-orbital 
contribution  to  the  various  energy  bands  of  NiS^  calculated  using  the 
interpolation  scheme  agrees  very  well  with  the  first  principles  calculation 
of  Bylander,  et  al.1^  Thus.  the  eigenfunctions  of  the  interpolation 
11am  il  toni  an  are  also  good  approximations  to  the  first  principles  results, 
end  perhaps  can  be  used  in  the  calculation  of  various  one-electron 
properties. 

The  broadening  of  the  Jd-bands  of  NiSij  with  respect  to  Ni  is  caused 
by  hybridization  effects,  since  the  A  parameters  of  Table  I  show  that  the 
direct  d-d  interations  are  much  smaller  in  the  compound  than  in  the  metal. 
The  larger  lattice  constant  of  the  intermetal  1 ic  compound  has  the  effect  of 
pushing  the  plane  wave  bands  downward  energetically  and  increasing  the 
density  of  plane  wave  states  below  Eg.  relative  to  fee  Ni.  The  d-bands  are 
degenerate  with  some  of  these  additional  plane  waves,  and  the  mixing  that 
results  produces  d~like  states  that  are  strongly  split  at  f  and  highly 
dispersive. 

A  similar  comparison  of  the  d-band  widths  can  be  made  between  the  Au 
intermetal  1 ic  compounds  and  Au  by  using  the  interpolation  scheme  to 
generate  a  nonr el ativ istic  Au  band  structure.  The  5d  band  width 
EC X5 >  — E(X^ )  for  Au  is  4.9  eV,  while  equivalent  band  widths  for  AuA^,  AuGa^ 
and  Auln^  are  1.05,  1.40,  0.56  eV,  respectivly.  In  this  case,  the  larger 


lattice  constant*  of  the  compounds  clearly  lead  to  narrower  d-band  widths 
in  the  AW*'  calculations.1*  This  is  because  the  Sd-bands  essentially  appear 


in  a  ba  nd 

gap 

region 

of  the 

pi  ane 

wave  states. 

and 

do  not  hybridize 

appr  eci ably 

with 

them. 

as  may 

seen 

by  looking 

at 

the  d-orbital 

contributions  to  the  various  bands  in  Table  II.  Comparing  the  A  parameters 
of  Au  and  its  compounds  show  that  Aj  and  A^ ,  which  directly  affect  the 
d-band  dispersion  and  thus  the  width.  are  very  much  smaller  for  the 
compounds  than  for  Au.  The  S  or  thogonl  ization  term  and  the  B  hybridization 
term  for  the  fluorite  compounds  are  also  significantly  less  than  those  of 
the  fee  metals.  However,  as  observed  in  the  case  of  NiSi^,  this  does  not 
necessarily  imply  that  the  d-bands  and  plane  waves  interact  weakly.  The  B 
values  for  the  Au  compounds  are  generally  smaller  than  those  for  NiSi^,  but 
the  S  terms  are  approximately  the  same  for  all  four  intermetal  1 ics.  The 
importance  of  these  terms  in  the  compounds  as  compared  to  the  fee  elements 
comes  from  the  higher  plane  wave  state  density  in  the  energy  region  near 
the  d-bands.  The  values  of  and  B£  for  the  Au  alloys  are  quite  different 

from  one  another,  unlike  case  of  elemental  Au.  This  is  reflected  in  the 
smaller  dispersion  of  e^  bands  than  that  of  t2(  bands  for  the  Au  alloys. 

The  spin-orbit  splitting  in  the  Sd  transi tion  metal s  is  quite  large, 
bat  performing  a  first  principles  band  structure  calculation  including 
spin-orbit  and  other  relativistic  effects  is  a  demanding  task.  However, 
once  a  set  of  parameters  for  the  nonrel  ativ  istic  interpolation  Hamiltonian 
have  been  determined,  it  is  relatively  straightforward  to  generate  a  band 
structure  including  spin-orbit  couplng  using  Eq.  6,  which  only  requires  one 
additional  parameter.  The  spin-orbit  parameter  (  may  either  be  taken 
directly  from  spectroscopic  determinations  of  atomic  energy  levels,  or  it 


16 


nay  be  determined  explicitly  for  the  material  of  interest  if  ARPES  data 

22 

:i\tt  which  show  the  d-band  splittings  at  the  [  point  of  the  BZ.  Since 
ARPLS  data  exist  for  both  AuGa^  and  Auln^,**  the  latter  approach  was  nsed 

here. 

Ihe  ARPLS  data  **  showed  three  clearly  resolved  d-bands  at  f,  which 
could  be  analyzed  in  terms  of  the  fee  tight-binding  Hamiltonian  including 
spin-orbit  coupling  to  yield  the  values  of  (  for  AuGa^  and  Auln^  shown  in 
Table  I.  However,  this  analysis  also  showed  that  the  binding  energies  of 
the  e  and  t~  levels  at  f  determined  by  the  A W  calculations**  were 
incorrect.  Therefore,  to  generate  the  band  structures  shown  in  Figs.  S  and 
6,  the  experimentally  determined  (  value  was  used  and  the  BO  and  A 
parameters  of  the  interpolation  Ham  il  toni  an  were  adjusted  to  force  the  band 
structure  to  agree  with  the  three  experimentally  determined  d-bands  at  [• 
All  the  other  parameters  used  in  generating  the  optimized  AuGa^  and  Auln^ 
band  structures  were  same  as  those  used  in  Figs.  2  and  3.  respectively. 

The  newly  determined  d-bands  were  less  strongly  bound  and  much  broader  than 
14 

the  AW  bands.  A  rather  fascinating  observation  is  that  the  relative 

splittings  of  the  Sd  bands  at  f  are  very  similar  for  Au,  AuGa^,  and 

Auln^,1*  but  the  overall  d-band  width  still  decreases  for  this  series  since 

the  dispersion  of  the  bands  decrease s  with  increasing  lattice  constant. 

13 

However,  for  AuGa^  and  Auln^,  the  d-band  width  determined  experimentally 
is  substantially  larger  than  the  AW  calculations  suggest. 

A  common  method  of  determining  the  quality  of  a  band  structure 
calculation  is  to  compare  the  energy  dependence  of  the  TDOS  to  spectral 
features  in  a  valence  band  XPS  spectrum.  Since  the  TDOS  curve  for  NiSi.  in 


11  £•  1  agrees  in  almost  every  detail  with  that  calculated  by  Tersoff  and 

7  3 

HaJtann/''  the  procedure  for  generating  the  1D0S  using  the  interpolation 
scheme  as  (h  scribed  in  Sec,  II  may  be  used  with  confidence.  In  making 
direct  comparisons  with  IPS  data,  however,  the  Gaussian  broadening  in  the 
TDOS  curves  should  be  larger  than  the  0.3  eV  used  in  Figs.  1-6  in  order  to 
account  for  lifetime  and  in  striaental  resolution  effects.  Using  a 
broadening  of  1.0  eV  to  generate  the  NiSij  ^S  produced  .  curve  ,ith  one 
broad  feature  that  agreed  very  well  with  reported  IPS  valence  band 
spectra.  H  J  Using  the  same  1.0  eV  broadening  for  the  nonrel  a  tiv  isti  c 
bands  of  the  Au  intermetal  1  ic  compounds  in  Figs.  2-4  produced  H)OS  plots 
with  a  single  peak  in  the  valence  band  region  in  Fig.  7,  which  were  in 
severe  disagreement  with  the  corresponding  IPS  spectra  of  van  Attekum,  et 
al.19  The  experimental  IPS  d-band  features  were  broader,  lower  In  binding 
clergy,  and  had  a  distinct  overlapping  double  peak  diape.  Ibwever,  dien 
the  TDOS  corresponding  to  AuGa2  and  Auln2  bands  in  Figg.  5  and  6> 
respectively,  and  the  adjusted  AuA12  j,anda  were  generated  wi  th  1.0  eV 
broadening  as  iiown  in  Fig.  7,  the  agreement  with  the  corresponding  IPS 
spectra  pas  very  good.  The  double  peak  is  apparent  in  the  adjusted  1DOS, 
and  the  relative  heights  of  the  two  peaks  agree  well  for  all  three 
compounds.  The  detailed  differences  between  the  IPS  da  ts  and  the  TDOS 
curves  may  be  the  result  of  photoelectron  cross  section  effects. 


V.  Conciliation 


The  sized  interpolation  scheme  for  fee  d-band  aetals  was  extended  to 
include  iluonte  structure  compounds  by  using  39  plane  waves  in  the  basis 
set.  The  band  structures  for  NiSi^,  AuAlj,  AuHa^,  and  Auln^  produced  by 
this  scheme  were  in  good  agreement  with  the  corresponding  first  principles 
calculations  up  to  3  eV  above  the  Feral  level,  with  ras  deviations  between 
the  first  principles  and  interpolated  bands  ~0.1  eV.  The  fit  to  the  higher 
energy  bands  was  not  as  good,  but  can  be  iaproved  by  using  more  plane  waves 
in  the  basis  set,  aore  terms  in  the  Haail toni an,  or  perhaps  a  better 
approximation  to  the  plane  wave  -  d  overlap  integrala.  Using  the 
interpolation  scheme,  spin-orbit  effects  were  added  in  a  very  direct  Banner 
to  a  nonrelativistic  band  structure  calculation.  Also,  it  was  possible  to 
adjust  the  parameters  of  the  interpolation  Haailtonian  to  obtain  agreeaent 
of  the  calculated  bands  at  f  of  the  BZ  with  AuGa^  and  Auln^  ARPES  data. 
The  TDOS  generated  using  the  interpolated  bands  of  NiSi^  and  the 
empirically  adjusted  bands  of  AuAl^,  AuGa^  and  Auln£  agreed  very  well  with 
XPS  valence  band  spectra.  In  the  case  of  the  Au  interne tal 1  ic  compounds, 
the  adjusted  d-bands  were  substantially  different  from  those  reported  by 
Switendick  and  Narath,^  which  did  not  agree  with  the  XPS  results. 
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Overlap  integrals  with  Slater  orbitals. 

The  radial  parts  of  d-wave  functions  RQ(j(r)  ar e  approximated  with 
;1  e  zeta  Slater  wave  functions, 


^nd ( r  ) 


cr 


n-le-?r. 


(Al) 


resulting  radial  overlap  integrals  are  represented  as  follows. 


^jj0-2  ( k  r )  t *  t~  frdr 
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Symmetrizing  factor* 

The  purpose  of  the  syaaetrizing  factors  is  to  compensate  for  the  use 
a  truncated  basis  set.  In  this  paper,  the  symmetrizing  factor  F. 


responding  to  plane  wave  k  has  the  following  form: 
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F 

1 


0 

l-(k2-560)2/19600. 


if  k^  <  560 
if  k2  >  700 

i 

otherwise. 


The  units  of  k^  are  such  that  the  I  point  in  the  first  B.Z.  is  at  (0.0,8), 
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ps  in  Smith's  scheme.  In  the  region  of  560<kt<700,  the  symmetrizing 
lectors  vary  smoothly  from  0  to  1 . 


Ref s. [ 14 ,1 5 , 17 ,20]  .  In  this  scheme  each 
used  as  an  offset  for  the  energy  levels. 


fable  1.  Parameters  used  in  the  sued  interpolation  icbati 


Table  11.  The  d-orbi tal  contribution  to  the  wave  function  of  the  bands  at  f  and  X 
I'.':  ct  s  (uni  t  s  in  %)  . 


-7.96 

14 

13.5 

-7.94 

91.9 

-7.52 

95.5 

-6.95 

99.1 

-3  .40 

74 

66.8 

-4.19 

6.8 

-3.80 

1.6 

-3.96 

0.8 

1.49 

32 

19.2 

1.08 

1.3 

2.17 

0.25 

1.11 

0.1 

(a)  energy  referenced  with  respect  to  Ep 

(b)  from  Ref. (IS] 


Figure  Captions 


(Fig.  1)  Energy  bands  of  NiSi^  generated  using  the  paraaieters  of  Table  I 
and  plotted  along  high  synnetry  lines  in  the  fee  Brillouin  Zone.  Also 
shown  is  the  corresponding  TDOS.  The  energy  scale  is  referenced  to  the 
Fermi  energy  and  TDOS  scale  has  arbitrary  units. 

(Fig.  2)  Energy  bands  and  TDOS  of  AuAl^  with  the  sane  renarks  as  Fig.  1. 

(Fig.  3)  Energy  bands  and  TDOS  of  AuGa^  with  the  sane  renarks  as  Fig.  1. 

(Fig.  4)  Energy  bands  and  TDOS  of  Auln^  with  the  sane  renarks  as  Fig.  1. 

(Fig.  5)  Energy  bands  and  TDOS  of  AuGa^  with  spin-orbit  interactions  in 

the  d-bands  and  adjnstnents  to  agree  with  the  ARPES  data  of  Ref.  13. 

(Fig.  6)  Energy  bands  and  TDOS  of  Aol^  with  the  aane  renarks  as  Fig.  3. 

(Fig.  7)  Conparison  between  the  TDOS  for  AuAl^,  AnGa^.  and  Auli^  end  the 
corresponding  XPS  spectra  of  Kef.  119],  all  plotted  versus  binding  energy 
referenced  to  the  Ferni  level.  XPS  spectra  are  shown  as  circles,  and  only 
represent  occupied  states.  TDOS  corresponding  to  nonrelativ istic  and 
relativistic  energy  bands  are  represented  as  dashed  and  solid  lines, 
respectively,  and  represent  both  occupied  and  enpty  states. 
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